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Abstract 

We calculate the evolution of quark-gluon-plasma droplets during the 
hadronization in a thermodynamical model. It is speculated that cooling 
as well as strangeness enrichment allow for the formation of strangelets even 
at very high initial entropy per baryon S'/^™'' ~ 500 and low initial baryon 
numbers of ^g'' ~ 30. It is shown that the droplet with vanishing initial 
chemical potential of strange quarks and a very moderate chemical potential 
of up/down quarks immediately charges up with strangeness. Baryon densi- 
ties of ~ 2po and strange chemical potentials of /Ug > 350 MeV are reached 
if strangelets are stable. The importance of net-baryon and net-strangeness 
fluctuations for the possible strangelet formation at RHIC and LHC is em- 
phasized. 
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Strangelets can be thought of as strange multiquark clusters which should be more com- 
pressed than ordinary nuclei and may exist as (meta-) stable exotic isomers of nuclear matter 
||1|]. It was speculated [0 that strange matter might resolve the dark matter issue, if it would 
be absolutely stable. 

The possible creation — in heavy ion collisions — of such long-lived remnants of the 
quark-gluon-plasma, cooled and charged up with strangeness by the emission of pions and 
kaons, was proposed by Liu and Shaw P] and Greiner et al. @,|^. Thus, strangelets can 
serve as unambiguous signatures for the creation of a quark gluon plasma. The detection 
of strangelets would verify exciting theoretical ideas with consequenses for our knowledge of 
the evolution of the early universe, the dynamics of supernova explosions and the underlying 
theory of strong interactions 0. 

Here we want to point out that such exotic states of matter can be created in heavy ion 
collisions even at collider energies, where such a process has received no attention so far, 
because common belief was that the (strange) baryon densities vanish at midrapidity, both 
at RHIC and LHC. We argue, however, that this conclusion was premature. This is due to 
the following effects: 

• fluctuations of the stopping power can provide finite baryochemical potential /ie at 
mid-rapidity in a small fraction of all events; 

• fluctuations of the net-baryon and -strangeness content between different rapidity bins 
within one event can be quite large; 

• strange (anti-)baryon enhancement due to collective effects (e. g. a chiral phase tran- 
sition); 

• strangeness and baryon distillery, which are inherent for the two-phase system (hadron 
gas/quark gluon plasma) for a wide parameter range. 

The last point stresses the significance of the 'chemistry' of the system in the evolution of 
the phase transition. 
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In the following we adopt a model for the dynamical creation of strangelets via the 
strangeness separation mechanism [Q. Consider a first order phase transition of the QGP 
to hadron gas. Strange and antistrange quarks do not hadronize at the same time for 
a baryon-rich system Q]. The separation mechanism can be viewed as being due to the 
associated production of kaons (containing s quarks) in the hadron phase, because of the 
surplus of massless quarks compared to their antiquarks. The strange quarks can combine 
to A-particles, but it is energetically favourable that s-quarks remain in the plasma, when 
hadronization proceeds. The ratio fs of the net strangess over the net baryon number 
quantifies the excess of net strangeness. Both the hadronic and the quark matter phases 
enter the strange sector 7^ of the phase diagram almost immediately, which has up to 
now been neglected in almost all calculations of the time evolution of the system. Earlier 
studies addressed the case of a baryon-rich QGP with rather moderate entropy per baryon 
IQJ^J^. Now we focus on low initial baryon densities and high specific entropies, to match 
the expected conditions of heavy ion collisions at RHIC and LHC, where the search for 
strangeness enters in the objective of the ALICE experiment [^. 

The hadronization transition has been described by geometric and statistical models, 
where the matter is assumed to be in partial or complete equilibrium during the whole 
( quasi- )isentropic expansion. A more realistic scenario must take into account the particle 
radiation from the surface of the hadronic fireball before 'freeze out'. Our model com- 
bines these two pictures. The expansion of the QGP droplet during the phase transition 
is described as a two-phase equilibrium; in particular the strangeness degree of freedom 
stays in chemical equilibrium because the complete hadronic particle production is driven 
by the plasma phase. The nonequilibrium radiation is incorporated by rapid freeze-out of 
hadrons from the outer layer of the hadron phase surrounding the QGP droplet. During the 
expansion, the volume increase of the system thus competes with the decrease due to the 
freeze-out. The global properties like (decreasing) S/A and (increasing) fs of the remaining 
two-phase system then change in time according to the following differential equations for 
the baryon number, the entropy, and the net strangeness number of the total system: 
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— = — r A^'^ 

dt 

^^*''* = -r5^« (1) 
^(iv,-iv^)*°* = -r(iv,-iv^)^^ , 

where T = -^j^ ( ^^t° ) effective ('universal') rate of particles (of converted liadron 

gas volume) evaporated from the hadron phase. The equation of state consists of the bag 
model for the quark gluon plasma and a mixture of relativistic Bose-Einstein and Fermi- 
Dirac gases of well established strange and non-strange hadrons up to in Hagedorn's eigen- 
volume correction for the hadron matter 0]. Thus, one solves simultaneously the equations 
of motion (|l]) and the Gibbs phase equilibrium conditions for the intrinsic variables, i.e. the 
chemical potentials and the temperature, as functions of time. 

Fig. illustrates the increase of haryon concentration in the plasma droplet as an in- 
herent feature of the dynamics of the phase transition (cf. [0]). The origin of this result 
lies in the fact that the baryon number in the quark-gluon phase is carried by quarks with 
mq <^ Tc, while the baryon density in the hadron phase is suppressed by a Boltzmann 
factor exp(— mbaryon/^c) with mbaryon ^ ^c- Mainly mesons (pions and kaons) are created 
in the hadronic phase. More relative entropy S/A than baryon number is carried away 
in the hadronization and evaporation process i.e. {S/A)^'-^ ^ [S / A)^^^ . Ultimately, 
whether {S/A)^^ is larger or smaller than [S/A)^^^ at finite, nonvanishing chemical po- 
tentials might theoretically only be proven rigorously by lattice gauge calculations in the 
furure. However, model equations of state do suggest such a behaviour, which would open 
such intriguing possibilities as baryon inhomogenities in ultrarelativistic heavy ion collisions. 
In the early universe shrinking quark droplets may — in analogy — contain the accumu- 
lated baryon number with possibly very high baryon density [Q]. This mechanism yields a 
primeval inhomegeneous nucleosynthesis which is signaled by the abundances of the 

light elements. 

What 'initial' conditions do we expect at collider energies? At RHIC energies one might 
see baryon stopping, dNB/dy > 0, on the average, at midrapidity. This can be due to 



multiple rescattering, leading to a nonvanishing, positive quarkchemical potential fig [0] . On 



the other hand, relativistic meson-field models, which, at high temperature, qualitatitively 
simulate chiral behaviour of the nuclear matter, exhibit a transition into a phase of massless 



baryons [|Tl|]. Including hyperons and yF-interaction |jT2[ it shows that at /i ^ the 
densities for all baryon species are of the order of po near the critical temperature. Thus, 
the fraction of (anti-) strange quarks increases drastically. Several hundred (anti-)baryons, 
many of them (anti-)hyperons may then fill the hot midrapidity region (with net baryon 
number ^ 0). Districts of non-vanishing net baryon (respectively anti-baryon) density with 
finite s (s) content will then occur stochastically. Thus, the finite chemical potential is locally 
caused by the fiuctuations of newly produced particles, not by the stopped matter. If such 
a phenomenon persists also for the deconfined phase, the effect of baryon concentration and 
strangeness separation may then result in the production of strangelets and anti-strangelets 
in roughly equal amounts. 

Let us try to give a rough estimate of the possible size of fluctuations in the net baryon 
and net strangeness number to be expected at midrapidity (or fluctuations along different 
rapidity intervals) around their mean. The average number of initial quarks and antiquarks 
(before hadronization) in a rapidity interval is approximately l/2dNT^/ dy ■ Ay, if half of the 
pions are made by the quarks and the other half by the gluons. For RHIC energies dNT^/dy 
has been estimated to lie between 1100 and 1600 and for LHC energies to lie up to nearly 



4000 in central Au+Au collisions ||T3|. Hence the quark number is roughly 500 for RHIC 



and up to 2000 for LHC in a rapidity interval Ay ~ 1. (In an equilibrated plasma the 
total number of quarks is ~ PgirR'^Az within Az = 1 — 2 fm in the early stage of the 
hydrodynamical expansion. According to pg = g^T^({3) ~ I.IT^ for a degeneracy of g=12 
these numbers correspond to temperatures of T ~ 250 — 500 MeV.) We take the number to 
be A^ = 500. A similar consideration holds for strange and antistrange quarks, and we take 
here Ns = 200. We now assume independent fluctuations according to Poissonians within 
this rapidity interval. In fact the actual width of the fluctuation at collider energies could be 
much broader (KNO-scaling of particle multiplicity distributions in elementary pp-coUisions 



P^). To justify the assumption of independent fluctuations of B and B despite of the local 
compensation of quantum numbers, one has to estimate the typical relative momenta within 
a quark-ant iquark pair. If one follows the parton cascade concept embodying perturbative 
QCD |T5[, the average of first parton-parton interactions {gg qq being the most 
important contribution) should be of the order of 5 — 10 GeV at LHC energies. The produced 
B and B, carrying about 0.4 of the (anti-)quark momenta, would thus be separated in 
rapidity by at least one unit (assuming a transverse momentum of about 500 MeV). 

The net baryon number in the box described above will be > 30 with a probability of 
0.5 %. About 0.1 % of the events will reach \B\ > 30 with a strangeness fraction |/s| > 0.7. 
Hence fiuctuations are not negligible. If each pion carries about 3.6 units of entropy (which 
is true for massless bosons), the entropy per baryon content in the fireball is 

As ^ dNs/dy ' 

and thus for dNs/dy = 30 a range of 60 to 250 is formed. If the plasma is equilibrated, 
the ratio of the quarkchemical potential and the temperature |/i|/T is directly related to the 
entropy per baryon number via 



o \ /I I \ — 1 

~ -7]- ' 



As\ 15 



(3) 



Accordingly the ratio then varies between 0.1 to 0.4. 

We now consider various fireballs with an initial net baryon number oi As = 30 and 
a net strangeness fraction fg of either or 0.7. The initial entropy per baryon ratios are 
chosen between 50 and 500. Table 1 summarizes the initial conditions (adjusted by the 
(initial) chemical potentials /ig, /i^ and temperature) used to start the hadronization. It also 
shows the final parameters of the quark droplet like the saturated strangeness content and 
baryon number. One further, yet crucial input, is the bag constant employed to describe 
the equation of state of the (strange) quark matter droplet. Only for the bag constants 
S^/^ < 180 MeV strange matter does exist as a metastable state at zero temperature [Q, 
being absolutely stable only for B^^^ < 150 MeV 13]. 



Fig. shows the time evolution of the baryon number for S/A^™^ = 200 and /™'* = 0.7 for 
varyous bag constants. For B^^^ < 180 MeV a cold strangelet emerges from the expansion 
and evaporation process, while the droplet completely hadronizes for bag constants B^^^ > 
180 MeV (for B^^^ = 210 MeV hadronization proceeds without any significant cooling of 
the quark phase, although the specific entropy S/A decreases by a factor of 2 from 200 to 
only 100). The strangeness separation works also in these cases, as can be read off the large 
final values of the net strangeness content, fs ~ 1.5 — 2. However, then the volume of the 
drop becomes small, it decays and the strange quarks hadronize into A-particles and other 
strange hadrons. 

Fig. 1^ shows the evolution of the two-phase system for S'/A'"* = 200, /™'* = and for a 
bag constant B^^^ = 160 MeV in the plane of the strangeness fraction vs. the baryon density. 
The baryon density increases by more than one order of magnitude! Correspondingly, the 
chemical potential rises as drastically during the evolution, namely from /i* = 16 MeV to 
jjf > 200 MeV. The stran geness separation mechanism drives the chemical potential of the 
strange quarks from /i* = up to /i{ ~ 400 MeV. Thus, the thermodynamical and chemical 
properties during the time evolution are quite different from the initial conditions of the 
system. 

Even for high initial entropies, S/A ^ 100 — 500, in the quark blob the entropy in the 
remaining droplet approaches zero at the end of the evolution (assuming B^/^ = 160 MeV). 
High initial entropies per baryon require more time for kaon and pion evaporation in order 
to end up finally with the same configuration of (meta-) stable strange quark matter. 

In conclusion, we have shown in the present model that the evolution of quark-gluon- 
plasma droplets during their hadronization may result in the formation of strangelets even 
at very high initial entropy per baryon 5/^4™'* ^ 500 and low initial baryon numbers of 
Ag'* ^ 30. The distillation of very small strangelets of a size As < 10 (see Table 1) is 
possible. We note that finite size effects of describing small strangelets neglected here might 
become of crucial importance ||T^. Special (meta-)stable candidates are the quark-alpha 



T7| with Ab = Q and the H-Dibaryon state with Ab = 2 |jT8[. Local net-baryon and net- 



strangeness fluctuations can provide suitable initial conditions for the possible strangelet 
creation at RHIC and LHC. Droplets with vanishing initial chemical potential of strange 
quarks and a small chemical potential of up/down quarks quickly charge up with strangeness 
and baryon-number: if strangelets are stable, the droplet reaches strange chemical potentials 
of Us > 350 MeV and two times ground state nuclear matter density! 



ACKNOWLEDGMENTS 



C.G. thanks the Alexander von Humboldt Stiftung for his support by a Feodor Lynen 
scholarship. This work was supported in part by the U.S. Department of Energy (grant DE- 
FG05-90ER40592), by the Gesellschaft fiir Schwerionenforschung, Darmstadt, Germany, the 
Bundesministerium fiir Forschung und Technologic (F. R. G.), Deutsche Forschungsgemein- 
schaft and the GraduiertenkoUeg Schwerionenphysik. 



8 



REFERENCES 

[1] A. R. Bodmer, Phys. Rev. D4, 1601 (1971); S. A. Chin and A. K. Kerman, Phys. Rev. 
Lett. 43, 1292 (1979); J. D. Bjorken and L. D. McLerran, Phys. Rev. D20, 2353 (1979) 

[2] E. Witten, Phys. Rev. D30, 272 (1984); E. Farhi and R. L. Jaffe, Phys. Rev. D30, 2379 

(1984) 

[3] H.-C. Liu and G.L. Shaw, Phys. Rev. D30, 1137 (1984) 

[4] C. Greiner, P. Koch and H. Stocker, Phys. Rev. Lett. 58, 1825 (1987); C. Greiner, D. H. 
Rischke, H. Stocker and P. Koch, Phys. Rev. D38, 2797 (1988) 

[5] C. Greiner and H. Stocker, Phys. Rev. D44, 3517 (1992) 

[6] Proc. Int. Workshop on Strange Quark Matter in Phys. and Astr., Nucl. Phys. (Proc. 
Suppl.) 24B, 1-299 (1991), see H. Kurki-Suonio 67, and S. Bonometto 54 

[7] K. S. Lee, U. Heinz, Phys. Rev. D47, 2068 (1993) 

[8] Letter of Intent, CERN/LHCC/93-16/LHCC/I 4; J. P. Coffin, project presentation, 
meeting of the ALICE technical board, April 23 (1995) 

[9] J. H. Applegate and C. J. Hogan, Phys. Rev. D31, 3037 (1985) 

[10] A. Dumitru, D. H. Rischke, T. Schonfeld, L. Winckelmann, H. Stocker and W. Greiner, 
Phys. Rev. Lett. 70, 2860 (1993) 

[11] J. Theis, G. Graebner, G. Buchwald, J. Maruhn, W. Greiner, H. Stocker, J. Polonyi, 
Phys. Rev. D28 (1983) 2286 

[12] J. Schaffner, C. Greiner, C. B. Dover, A. Gal, H. Stocker, Phys. Rev. Lett. 71 (1993) 
1328 

[13] X.-N. Wang and M. Gyulassy, Phys. Rev. D.44, 3501 (1991); K. Geiger, Phys. Rev. 
D47, 133 (1992) 



9 



[14] UA5 coUab., G. J. Alner et al., Phys. Rep. 154 (1987) 247 
[15] K. Geiger and B. Miiller, Nucl. Phys. B 369 (1992) 600 

[16] J. Madsen, Phys. Rev. Lett. 70, 391 (1993) and Conf. on Strangeness and Quark Matter, 
Crete, Greece (1994) 

[17] F. C. Michel, Phys. Rev. Lett. 60, 677 (1988) 

[18] R. L. Jaffe, Phys. Rev. Lett. 38, 195 (1977) 



10 



FIGURES 




FIG. 1. Time evolution of the net baryon density of a QGP droplet. The initial conditions are 
/]°" = and A^'^ = 30. The bag constant is = leo MeV. 
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TABLES 

TABLE I. Various situations of an hadronizing plasma droplet. In the first column the bag 

constant for describing the plasma phase is listed. Then the initial conditions follow. The final 
values for the baryon number, the strangeness fraction and the two chemical potentials at the end 
(or after) hadronization are listed in the last four columns. 
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